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Nuclear physics is a vast and rich field whose phe- 
nomenology has been explored for decades through in- 
tense experimental and theoretical effort. However, there 
is still little understanding of the connection to the under- 
lying theory of the strong interactions, Quantum Chro- 
modynamics (QCD), and the basic building blocks of 
nature, quarks and gluons. A key element in connecting 
QCD to nuclei is determining how the properties of the 
hadrons, and their responses to external probes, can be 
described and predicted in terms of quark and gluon dis- 
tributions. This is the focus of both theoretical and ex- 
perimental efforts world-wide, and in particular, at the 
Jefferson Laboratory. At the next length scale up from 
single hadrons are nuclei, which are systems of neutrons 
and protons (and also hyperons) that are bound by the nu- 
clear forces, and it is this area that traditionally defined 
the field of nuclear physics. Nuclear physics also encom- 
passes the study of matter under extreme conditions of 
temperature and density. This continues to be the subject 
of extensive experimental efforts using heavy-ion colli- 
sions, such as produced at RHIC at Brookhaven National 
Laboratory and now at the LHC at CERN. The range of 
densities that are explored in such collisions is not great, 
and does not match those found in some stellar environ- 
ments such as in the late-stages of evolution of Type II 
Supernova - the birthplace of the heavy-elements. Nu- 
clear physics also impacts our understanding of funda- 
mental cosmological aspects of our universe. Type I Su- 
pernova are used as standard candles to 'calibrate' the 
universe. In fact, they provided the first measurement of 
the acceleration of the universe and presently provide a 
good constraint on the dark-energy-dark-matter compo- 
sition of the universe. 

During 2009, a series of workshops were held in the 
United States that established the scientific needs for 
exa-scale computational resources for the field of Nu- 
clear Physics. Just to set the scale - 1 Exaflop = 1000 



Petaflops =10^ Teraflops = 10^ Gigaflops. Partner work- 
shops were held in the areas of Biology, Basic Energy 
Sciences, Climate Science, Fusion Energy, High Energy 
Physics, National Security, Nuclear Energy, Architec- 
tures and Technology, and Cross-cutting areas, the details 
of which can be found in Ref. [1]. In 2010 the output, 
conclusions, and recommendations from the panels at the 
workshops were synthesized into an overall report. The 
nuclear physics panel determined that exa-scale comput- 
ing is required for the calculation, with quantifiable un- 
certainties, of systems and processes involving strongly 
interacting matter over a large range of temperatures and 
densities. 

Due to the increasing computational resources becom- 
ing available to the field, progress is being made toward 
a direct connection between QCD and nuclear physics 
using the numerical technique of Lattice QCD (LQCD). 
It has been known for the past four decades that QCD, 
together with the electroweak interactions, underlies all 
of nuclear physics. However, soon after the discovery of 
QCD it became apparent that the complexities of the the- 
ory at strong coupling would hinder analytic progress in 
understanding the properties of the simplest hadrons, let 
alone the simplest features of the nuclear forces. Wil- 
son pointed the way to eventual direct quantitative con- 
firmation of the origins of nuclear physics by formulat- 
ing Lattice QCD (LQCD) [2], a regularization and non- 
perturbative definition of QCD that is suitable for the in- 
tensive computational demands of solving QCD in the 
infrared. LQCD is a technique in which Euclidean space 
correlation functions are calculated by a Monte-Carlo 
evaluation of the Euclidean space path integral [2]. The 
calculations are performed in Euclidean space so that 
field configurations that have a large action are expo- 
nentially suppressed. This is in contrast with Minkowski 
space in which large action contributions result in a com- 
plex phase which will average to an exponentially small 



contribution with nearby configurations. Space-time is 
discretized with the quarks residing on the lattice sites, 
and the gluon fields residing on the links between lat- 
tice sites. The lattice spacing, the distance between ad- 
jacent lattice sites, is required to be much smaller than 
the characteristic hadronic length scale of the system 
under study. The effects of a finite lattice spacing can 
be systematically removed by combining calculations 
of correlation functions at several lattice spacings with 
the low-energy effective field theory (EFT) which ex- 
plicitly includes the discretization effects. This type of 
EFT is somewhat more complicated than its continuum 
counterpart as it must reproduce matrix elements of the 
Symanzik action constructed with higher dimension op- 
erators induced by the discretization [3]. While the ac- 
tion lacks Lorentz invariance and rotational symmetry, 
it is constrained by hypercubic symmetry. As comput- 
ers have finite memory and performance, the lattice vol- 
umes are finite in all four space-time directions. Gener- 
ally, periodic boundary conditions (EC's) are imposed 
on the fields in the space-directions (a three-dimensional 
torus), while (anti-) periodic EC's are imposed on the 
(quark) gauge fields in the time-direction, which in many 
cases is much larger than the space-directions (in or- 
der to approach the zero-temperature limit). During the 
last five years, LQCD has emerged from a long period 
of research and development to the point where pre- 
cise predictions for hadronic quantities are starting to be 
made. In particular, presently, fully-dynamical calcula- 
tions with near-exact chiral symmetry at finite lattice- 
spacing have become standard, with lattice volumes of 
spatial extent 2.5 fm and with lattice spacings in the 
range 0.12 fm. It is still the norm that the light-quark 
masses, m^, are larger than those of nature, with typi- 
cal pion masses rriji^ 300 MeV. However, exploratory 
calculations at, and near, the physical pion mass, rriji ~ 
139 MeV, have been performed, but the lattice volumes 
at these quark masses remain too small for the results to 
be meaningfully compared with nature. 

Zero-temperature LQCD is a very active field. This 
area is focused on the determination of the properties of 
the hadrons, establishing the ground-state masses and ex- 
citation spectrum, and mapping out the parton and gener- 
alized parton distributions (as an example, see figure 1). 
This program is closely tied to the parallel experimen- 
tal efforts, in particular at the Jefferson laboratory, and 
has the exploration of exotic hadrons as a priority. Exotic 
hadrons are ones that are not predicted by the most naive 
quark-model, but follow naturally from QCD - states 
with interpolating operators that can be constructed only 
by including the gluon field. A main focus of the efforts 
in fundamental symmetries is the prediction of the neu- 
tron electric-dipole moment resulting from a non-zero 
value of 6, and from higher dimension operators that 



may also contribute. An estimate of the computational 
resources required to calculate the structure of the nu- 
cleon can be seen in figure 2. 
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FIGURE 1. LQCD calculadons of the first moment of the 
isovector parton distribution in the proton performed with the 
domain- wall fermion discretization of the quark-fields [4]. I 
thank the authors for allowing me to reproduce this figure. 
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FIGURE 2. Estimates of the computational resources re- 
quired to calculate the structure of nucleons [5]. 

When it comes to calculating the properties of the 
light nuclei and their interactions with LQCD, the first 
"benchmarking" step is for LQCD to post-diet the large 
nucleon-nucleon scattering lengths and the existence of 
the deuteron, the simplest nucleus. This will firmly estab- 
lish the connection between QCD and nuclear physics, 
and will allow for an exploration of how nuclei and 
nuclear interactions depend upon the fundamental pa- 
rameters of nature. In particular, it is believed that an 
understanding of the fine-tunings that permeate nuclear 
physics will finally be translated into fine-tunings of the 
light-quark masses. While these issues are of great in- 
terest, and it is important to recover what is known ex- 
perimentally to high precision, these goals are not the 



main objective of the LQCD effort in nuclear physics. 
The primary reason for investing resources in this area 
is to be able to calculate physical quantities of impor- 
tance that cannot be accessed experimentally, or which 
can be measured with only limited precision in the lab- 
oratory. Two important examples of how LQCD calcula- 
tions can impact nuclear physics are in the structure of 
nuclei and in the behavior of hadronic matter at densi- 
ties beyond that of nuclear matter. LQCD will be able to 
calculate the interactions of multiple nucleons, bound or 
unbound, in the same way that it can be used to determine 
the two-body scattering parameters. While it is not pos- 
sible to extract objects from LQCD calculations that can 
be directly compared with the modern nucleon-nucleon 
potentials (despite claims to the contrary), the scattering- 
phase shifts can be determined in the same way that 
experimentally measured cross-sections are translated 
into phase-shifts. LQCD calculations of three-nucleon 
systems will allow for the determination of the three- 
nucleon interactions. Steps towards this goal have been 
taken with analogous calculations in meson-systems, e.g. 
Ref. [6], and in the first calculations of three- and four- 
nucleon systems [7, 8]. Two decades of effort in the de- 
velopment of low-energy effective field theories (EFT's) 
for nuclear physics provide a framework with which to 
classify the contributions from multi-nucleon operators, 
and LQCD calculations will first match onto such oper- 
ators. This will provide some sort of unification for nu- 
clear physics as the tools will then be in place, includ- 
ing those used to perform nuclear structure calculations 
starting from modern nucleon-nucleon potentials, to start 
from QCD and calculate the properties of multi-nucleon 
systems, a long-standing goal of the field. The results 
of a recent calculation of baryon-baryon phase-shifts [9] 
by the NPLQCD collaboration can be seen in figure 3, 
and the estimated resource requirements to provide pre- 
cise calculations at the physical pion mass are shown in 
figure 4. Simulations of LQCD calculations have been 
performed [10] to determine the level of precision re- 
quired to extract the deuteron binding energy, the results 
of which can be seen in figures 5 and 6. 



In addition to the large computational resources and 
algorithmic improvements that are required to complete 
the mission, formal developments are also required. For 
instance, a reliable method with which to extract inelastic 
scattering cross- sections from LQCD calculations does 
not yet exist, and must be developed. Further, a bet- 
ter understanding of the convergence of EFT's that are 
used in such systems is required. Current LQCD calcu- 
lations indicate unexpected behaviors in the convergence 
patterns of the EFT's, including chiral perturbation the- 
ory (;t^PT). Two striking examples of such behaviors are 
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FIGURE 3. The most recent calculation [9] of baryon- 
baryon scattering phase-shifts at a pion mass of mji ~ 
390 MeV, calculated on lattices generated by the Hadron Spec- 
trum Collaboration [11, 12] with a spatial-lattice-spacing of 
a ^ 0.123 fm, lattice-volume of 20^ x 128 with an anisotropy 
factor of = 3.5. 
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FIGURE 4. Estimates of the computational resources re- 
quired to calculate the interactions between hadrons [5]. 




FIGURE 5. Nucleon-nucleon scattering parameters recon- 
structed from simulated LQCD calculations [10]. The left panel 
corresponds to 10% precision, the middle panel to 5% precision 
and the right panel to 1 % precision calculations of the energy- 
splitting between two interacting nucleons and two free nucle- 
ons. 



evident in LQCD calculations of the nucleon mass, ex- 
hibiting a light-quark mass dependence that is consistent 
with being linear in the pion mass (i.e. [13], as 

shown in figure 7, and in the meson-meson scattering 
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FIGURE 6. The deuteron binding energy reconstructed from 
simulated LQCD calculations [10]. The left panel corresponds 
to 10% precision, the middle panel to 5% precision and the 
right panel to 1% precision calculations of the energy- splitting 
between two interacting nucleons and two free nucleons. 



lengths which are consistent with their tree-level predic- 
tions, even at large light-quark masses [14]. 



lated into fine-tunings of the fundamental parameters. Of 
course, this requires calculating the structure of alpha- 
cluster states in ^Be and ^^C, which are notoriously diffi- 
cult starting from single particle-orbit type calculations. 
Such calculations are a major goal of the nuclear struc- 
ture community for the next decade, and preliminary cal- 
culations are underway. Nuclear structure calculations 
play an important role in determining the properties and 
interactions of neutrinos from present and future neutrino 
experiments. An estimate of the computational resources 
required to perform the nuclear structure calculations re- 
quired to analyze the results of neutrino experiments are 
shown in figure 9. 
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FIGURE 7. A summary of LQCD calculations of the mass 
of the nucleon) [13]. I thank A. Walker-Loud for allowing me 
to reproduce this figure. 

While the nuclear structure community has developed 
impressive tools to perform calculations starting from 
either effective interactions in truncated model- spaces, 
or in fully no-core- shell-model (NCSIM), large com- 
putational resources are required for such calculations. 
Like QCD, this involves quantum many-body systems, 
and exa- scale resources are required to perform calcula- 
tions of nuclei starting from the two-nucleon and multi- 
nucleon forces beyond the lightest few nuclei. An exam- 
ple of the state-of-the-art in the calculation of light nuclei 
starting from the nucleon-nucleon and the three-nucleon 
interaction [15] can be seen in figure 8 in which the role 
played by three-body nuclear forces is made clear. In ad- 
dition to being extremely important for predictive capa- 
bilities for nuclear processes in environments where ex- 
periments are impossible, such calculations will define 
how nuclei and their interactions depend upon the fun- 
damental parameters of nature, namely the strong inter- 
action length-scale, Aqcd, the quark masses, m^, and the 
electroweak interactions. The fine-tuning(s) required to 
produce sufficient carbon for us to exist will be trans- 



FIGURE 8. The spectrum of A = 6 - 12 nuclei calculated 
with GFMC [15]. I thank the authors for allowing me to repro- 
duce this figure. 
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FIGURE 9. Estimates of the computational resources re- 
quired to calculate the structure of nuclei required to extract 
neutrino properties from experimental measurements. 

The nature of calculations in heavy nuclei is quite dif- 
ferent to those in lighter nuclei. It is not presently fea- 
sible to start from the measured or calculated nucleon- 
nucleon interaction and the three-nucleon interaction and 



calculate, for instance, fission cross-sections directly. 
Calculations starting from an energy density functional 
have been a focus in the community for a number of 
years, e.g. Ref. [16]. Hartree-Fock-Bogoliubov calcula- 
tions [17] will require exa-scale computing to map out 
the energy- surface associated with the fission of heavy 
nuclei. An example of which is shown in figure 10, where 
it is clear that both cold and hot fission can occur. A 
multi-dimensional grid in the parameters describing the 
nuclear shape, Q20, Q22, Q30, Q40, and scission of the 
nucleus Qn, requires ~ 3 x 10^ points, and requires an 
exa-flop to execute one calculation in less than one day 
(see Ref. [18] for a more complete discussion). 



of the computational resources required to map out the 
QCD phase-diagram at zero-baryon-number density are 
shown in figure 13. 
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FIGURE 11. A cartoon of the evolution of a heavy-ion col- 
lision. 
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FIGURE 10. An example of an energy-surface that is 
found in the fission of a heavy-nucleus as function of shape- 
deformations 220 and 240- 1 thank W. Younes for allowing me 
to reproduce this figure. 

Exa-scale computing resources are required in order 
to determine the behavior of matter under extreme con- 
ditions of temperature and density. Enormous effort has 
already been placed into the exploration of equilibrium 
properties of hot matter, and to determining the behavior 
of matter during the collisions of relativistic heavy-ions 
(a cartoon of which is shown in figure 11). LQCD is the 
tool with which to calculate high-temperature properties 
of systems in equilibrium (see, for example, Ref. [19]), 
while the evolution of heavy-ion collisions requires a 
host of techniques, some of which are not yet in place. 
LQCD has established that the transition from confin- 
ing to deconfining at vanishing baryon number den- 
sity is a cross-over at the physical values of the light- 
quark masses, as opposed to, say, a first-order transi- 
tion, and significant effort is being placed on mapping 
out the transition in a number of observables with differ- 
ent quark discretizations, an example of which [20, 21] 
is the Polyakov-loop, as shown in figure 12. A signifi- 
cant conceptual problem that remains in this area is the 
inclusion of moderate and large baryon number density 
due to the infamous sign-problem that is encountered in 
numerical sampling of the partition function. Estimates 
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FIGURE 12. The Polyakov loop as a function of tempera- 
ture calculated with LQCD with a variety of quark discretiza- 
tions [20]. I thank the authors for allowing me to reproduce this 
figure. 
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FIGURE 13. Estimates of the computational resources re- 
quired to calculate the QCD phase-diagram at zero baryon 
number density. 

The other areas of nuclear physics that have been 
identified as requiring exa-scale computing resources to 



accomplish their goals are astrophysics and accelerator 
design [5]. These areas are beyond the scope of this 
meeting and I will not discuss them further. 

Nuclear physicists will need enhanced collabora- 
tion with applied Mathematicians, Computer Scientists, 
Statisticians, and others in order to prepare for, and to 
fully utilize, exa-scale computing resources when they 
become available - which is anticipated to occur around 
2017. Feedback has already been provided to those de- 
signing the machines by the community, and this will 
have to be an ongoing dialogue. As the nature of the ma- 
chines will be quite different from those of today, the 
codes and algorithms will need to evolve accordingly. 
The codes will have to become fault-tolerant - they will 
need to be able to survive the loss of multiple compute- 
cores. Further, as the memory-per-core will likely be sub- 
stantially less than in today's machines, the mode of cal- 
culation will have to change. The movement of informa- 
tion is a major consumer of power and will need to be 
minimized. Consequently, much of the post-processing, 
that is done today on smaller machines, will have to be 
done in situ. In many fields, visualization is crucial to ex- 
ploring the massive outputs from calculations. Nuclear 
and particle physics is somewhat lagging in this area - 
a situation that will have to change within the next few 
years. Given the scale of the computational facilities, one 
should view these numerical efforts like experiments in 
their style of operation. As the nuclear physics commu- 
nity could not efficiently use an exa-scale machine to- 
day, but could efficiently use a 1-10 sustained peta-flop 
machine, a staged evolution to the exa-scale will likely 
be required. However, it must be stressed that computing 
power is insufficient by itself to accomplish all of the 
goals. Conceptual and formal developments must take 
place in some areas, hand-in-hand, with the deployment 
of computing resources and developments in algorithms 
and applied mathematics. 

The field of nuclear physics will be somewhat rev- 
olutionized by the deployment of exa-scale comput- 
ing resources. Exa-scale computing will provide nuclear 
physics with predictive capabilities that allow for the as- 
signment of reliable uncertainty estimates in observables 
that cannot be explored experimentally. Further, it will 
enable the systematic exploration of fundamental aspects 
of nature that are manifested in the structure and interac- 
tions of nuclei. 
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